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ABSTRACT 

We explore the rate and impact of galaxy mergers on the massive galaxy population using the amplitude of the 
two-point correlation function on small scales for M, >5x 1O 1O M galaxies from the COSMOS and COMBO- 
17 surveys. Using a pair fraction derived from the Sloan Digital Sky Survey as a low-redshift benchmark, the 
large survey area at intermediate redshifts allows us to determine the evolution of the close pair fraction with 
unprecedented accuracy for a mass-selected sample: we find that the fraction of galaxies more massive than 
5 x 10 10 M Q in pairs separated by less than 30kpc in 3D space evolves as F(z) = (0.0130±0.0019) x (l+z) ul±a2S 
between z = and z= 1.2. Assuming a merger time scale of 0.5 Gyrs, the inferred merger rate is such that 
galaxies with mass in excess of 10 h Mq have undergone, on average, 0.5 (0.7) mergers involving progenitor 
galaxies both more massive than 5 x 10 10 M Q since z = 0.6 (1.2). We also study the number density evolution of 
massive red sequence galaxies using published luminosity functions and constraints on the M/Lg evolution from 
the fundamental plane. Moreover, we demonstrate that the measured merger rate of massive galaxies is sufficient 
to explain this observed number density evolution in massive red sequence galaxies since z = 1 . 

Subject headings: galaxies: general — galaxies: evolution — galaxies: interactions — galaxies: ellipticals and 
lenticulars, cD — galaxies: formation — galaxies: statistics 



1. INTRODUCTION 

The distribution of galaxies in a colour-magnitude diagram 
shows a bimodality: galaxies with red optical colors occupy a 
relatively tight sequence in color (the red sequence) while blue 
galaxies sho w a wider dispersio n and populate the so-called 
'blue cloud' dBlanton et alj|2003l) . This bimodal ity has existed 
for most of the history of the uni verse (e.g., [Bell et al. l2004t 
iFaber et al.l2007HTavlor et al.l2009t) . However, the stellar mass 
density of red sequence galaxies has increased by a factor ~ 2 
since z = 1, w hile it has remained approximately constant for 
blue g alaxies dBellet al J 120041: IFaber et al J 120071: iBrown et all 
2007 1). Yet, the ma j ority of stars are fo rmed in blue galax- 
ies dBell et all 120071: IWalcher et al.ll2008l) . which implies that 
galaxies migrate from the bl ue cloud to the red sequence (see, 
e.g., the discussion by IFaber et aT]|2007l) . 

While significant evolution of the integrated mass function of 
red galaxies is generally agreed upon, the evolution in the num- 
ber density of high- mass galaxies has been shown to be consis- 
tent with zero (e.g JCimatti et al]200"6l:ICool et al]2008l) . How- 
ever, given the large uncertainties, the constraints are not tight. 
There is ample room for considerable evolution, and such evo- 
lution is expected in a ACDM hierarchichal cosmology, which 
predicts significant growth of galaxies through mergi ng, espe- 
cially at the high-mass end dDe Lucia & Blaizotll2007l) . 

The most massive, non-starforming galax ies are essentially 
all sp heroidal (i.e., have ellipticities < 0.4; Ivan der Wel et al 



in their halos. 

The goal of this paper is to estimate the impact of galaxy 
mergers on the evolution of massive galaxies from z ~ 1 to the 
present day, and more specifically, on the massive end of the 
red sequence. To that end we measure the merger fraction since 
z = 1.2 in a large sample of mass-selected galaxies. We will 
test whether these measurements are consistent with predic- 
tions from hierarchichal formation models, and, providing an 
independent estimate of the number density of the most mas- 
sive red galaxies, with the observed evolution of the bright tail 
of the luminosity function of red galaxies. 

It is observationally challenging to identify galaxy merg- 
ers, especially at large cosmological distances. Mergers are 
found either in an early phase of the interaction, when the 
two galaxies have not yet coalesced and are found in a close 
pair dPatton et al ll2000[ iLe Fevre et ail 120001: iLin et al.ll2004 
2008i iBell et alj|2006al: kartaltepe et al.1 120071: IRobaina et all 



2009), or in a later phase, when they display signatures of 
gravita tional interaction and are just prior to or a f ter coales- 
cence d Abraham et alj|l996t IConselice et al.| 12003b lLotz et al.1 
2004 t IBell et alj|2005t Mcintosh et alJ|2008tj jogee et alj|2009h 
iHeiderman et al.1 120091 IRobaina et al] l2009h . A wide range 
of results have been found for the merger fraction evolution, 

parameterized as (1 +z) m , with m ranging from to 4 (e.g., 

Patton et aLlioOOtlLe Fevre et aL^OOOtLin et al.l2004l:lKartaltepe eTal] 



2007). 



onnm ■ a- , u ».u e a t u u ■ It- e t> "1 Here, we use robust 2-point correlation function techniques 

2009), indicating that they formed through merging (Toomre & Toomre ' . , . K , , _ . , „ M 

THRIFT — : — o o -4. Jinn-i LH^ j — i i l™rvnk ov. dti a sample of galaxies with 0.2 < z < 1.2 with stellar masses 



1972s lSchweizer & Seitzerll992l:lKormendv et al.l2009t) . Shape 
distributions do not specify when merging happened. Some 
valuable constraints have been placed by the study of the clus- 
tering properties of red galaxies in this mass range by Whit e et al] 
(2007), who found that 1/3 of massive red satellites disappear 
between z=0.9 and z=0.5 by merging with the central galaxies 



sample or galaxies witn u./ < z 
in excess of 5 x 1O 1O M selected from the COSMOS and 
COMBO-17 surveys to quantify the merger rate of massive 
galaxies and its evolution. We augment the statistical signifi- 
cance of our analysis by using an estimate of the pair fraction 
found in Sloan Digital Sky Survey (SDSS) at z ~ 0. 1 . The total 
volume probed by this study at intermediate redshifts is at least 
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4 times larger than any previous mass-selected study on the 
evolution of the merger fraction, and reduces dramatically the 
systematic uncertainties related to cosmic variance by the use of 
four independent fields. Then we compare the inferred galaxy 
merger rate with the observed number density evolution of mas- 
sive (M* > 10 n M Q ), red galaxies from z ~ 1 t o the present 
day, which we obtain by converting B rown et al. ld2007|)LFs to 
stellar-mass functions. We assume Hq = 70 km/s, Q, m = 0.3 and 
Q A = 0.7. 

2. SAMPLE AND METHOD 

The bulk of ou r sample is drawn f rom the ~2sq. degree 
COSMOS survey (IScoville et al.ll2007l) . We use photometric 
redsh ifts calculated from 30-band photometry by lllbert et al] 
(2009); comparison with spectroscopic redshifts shows excel- 
lent accuracy. We use those redshifts to derive rest-frame quan- 
tities and stellar masses by using the observed broad-band pho- 
tometry in conjunction with a non-evolving tem plate library de- 

ef - 



ll983HLi et al.l2006l) . the real-space 2pcf can be reasonably well 
fit by a power-law. Assuming £(r) = (r/Vo)~ 7 , then w(r p ) = 
Cr 7 ^" 7 , with C0F{r[( 7 - l)/2]/r( 7 /2)}. We fit the latter ex- 
pression to our data to find the parameters 7 and ro and use 
them in the real-space 2pcf to find the fraction of galaxies in 
close pairs. 

As we wish to preserve S/N we did not integrate along the 
entire l.o.s. when calculating vv(rp). Instead we allowed galax- 
ies to form a pair only if the redshift difference was smaller than 
Ppair = \/2 x cr., with c, being the redshift error of the primary 
galaxy. As the photo-z errors follow a Gaussian distribution 
with width a z dWolf et all 120031, |2004|) a fraction of the pairs 
would be missing by simply imposing the l.o.s. criteria. Thus, 
following lBeil et all (l2006al) the fraction of galaxy pairs is 



/ = 



2 /2c 



' dz. 



(2) 



rived using Pegase stellar po pulation model fsee lFioc & Rocca-Volmerarffifr i, w(rp) is multiplied by 1 // in order to account for miss- 



19991) and a lChabrierl (120 03) stellar initial mass function ( IMF). 



The use of a iKroupa et al.1 (1 19931) or a iKroupal (1200 lb IMF 
would yield similar stellar masses to within ~ 10%. The 
reddest templates are produced through single exponentially- 
declining star formation episodes, intermediate templates also 
have a low-level constant star formation rate (SFR) and the 
bluer templates have superimposed a recent burst of star forma- 
tion. A full description of the stel lar masses will be pr ovided in 
a future paper; comparison with IPannella et ai] (2009) masses 
shows agreement at the 0.1 dex level. 

To combat the sample variance of a single 2 sq. degree field, 
we augment the COSMOS sample with a sample drawn from 
three widely-separated 0.25 sq. degree fields from COMBO- 
17. COMBO-17 photo-z 's, colors and stellar masses have bee n 
extensively describ ed i n IWolf et all (l2003l)IWolf et all (|2004|) . 
iBorch et al l (120061) and lGrav et aTT~(12009l) . Given the different 
depths of the two surveys we only include galaxies from the 
COMBO-17 catalog at z < 0.8, where it is complete for our 
mass limit. The final sample comprises ~ 18000 galaxies with 
M, >5x 10 I0 M Q over an area of ~ 2.75 sq. deg. in the redshift 
range 0.2 < z < 1.2. 

We use the fraction of galaxies with a companion closer than 
30 kpc (close pairs) as a proxy for the merger fraction, as those 
systems are very likely to merge in a few hundred Myfl As 
redshift errors translate into line-of-sight (l.o.s.) distance un- 
certainties of the order of ^ 50-100 Mpc we use projected 
2-point correlation functions (2pcf) to find the number of pro- 
jected close pairs and then deproject into the 3D space. 

The projected correlation function w(rp) is the integral along 
the line of sight of the real-space correlation function: 



ing pairs; in our case a correction factor of 1 . 19. 

Given the 3D correlation function £(r), the differential prob- 
ability of finding a galaxy occupying a volume SV at a distance 
r of another galaxy is SP = n[l +£(r)]8V, where n is the number 
density of secondary galaxies. Then, by a simple integration 
of this expression, we obtain the probabili ty of a galaxy be- 
ing within a distance rf of any other g alaxy dPatton et alJ uOOO; 
iBell et al.ll2006atlMasiedi etai]|2006l) : 



Because £(r) 
obtain: 



w(r p )= / air 2 p + 7T 2 ]' /2 )d7r, 



(1) 



where r p is the distance between the two galaxies projected on 
the plane of sky and it the line-of-sight separation. A conve- 
nient and simple esti mator of the 2pcf at small scales is w(r D ) = 
A(DD/RR - 1) (e.g. IBell et all l2006at iRobaina et al] 120091) . 
where A is the path length being integrated over, DD(rp) is 
the histogram of separations between real galaxies and RR(rp) 
is the histogram of separations between galaxies i n a randomly- 
distributed catalog. As shown in the literature (e.g. Davis & Peebles 

'As the lower mass limit of our sample is 5 X 1O 1O M0, we are automatically 
selecting both members of the pair to be above that mass. 



P(r<r f )= / n[\+i(r)]dV (3) 
Jo 

= (r/ro)" 7 and £(r) ^ 1 at the small scales, we 



P(r <r/) = fpair = 



4™ 



"7 
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7 > 



(4) 



where / pfl ,> is the fraction of galaxies in close pairs. As galaxy 
interactions are completely decoupled from the Hubble flow, 
we calculate probabilities as a function of the proper (physical) 
separation between the two galaxies. Errors in the correlation 
function are calculated by means of bootstrap resampling. 

3. RESULTS AND DISCUSSION 

This work aims to an understanding of the impact of galaxy 
mergers on the evolution of massive galaxies and the evolution 
of the red sequence at its massive end. To that end we apply the 
above method to our sample of galaxies with > 5 x 1O 1O M , 
regardless of their optical colors, to find the fraction of galaxies 
in close physical pairs, which we will use as a proxy for the 
fraction of galaxies undergoing a galaxy interaction. Hereafter, 
all reference to close pairs will refer to galaxy pairs with in- 
ferred separations smaller than 30 kpc in 3D space. Fig. Q] and 
Table Q] show the result of our 2pcf analysis. 

The projected correlation functions in Fig.Q]are shown over 
the range 10 < r p /kpc < 1000, in physical (not comoving) 
units, as is appropriate for a process such as galaxy merging 
that has decoupled from the Hubble flow. One can see that the 
fitted power laws are good descriptions of the run of w{r p ) with 
radius, down to the resolution limit of ~ 15 kpc. Below this 
radius, we extrapolate the power law fit. We have checked the 
accuracy of this extrapolation by comparison with extremely 
close pairs resolved by HST in the s t udy of two of the four fields 
used in this paper bv lRobaina et al.l d2009h . An extrapolation of 
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Table 1 

3d correlation function parameters and close pair fractions 



z 


7 


ro 


n(Mpc J ) 


fpatr (< 30 kpc) 


0.2 < z < 0.4° 


2.03 ±0.05 


3.50±0.50 


0.0031 


0.0171 ±0.0050 


0.4 < z < 0.6" 


2.06 ±0.04 


3.60 ±0.50 


0.0034 


0.0238 ±0.0043 


0.6 < z < 0.8" 


1.94 ±0.04 


3.80±0.30 


0.0064 


0.0241 ±0.0038 


0.8 <z< 1.0* 


1.93 ±0.03 


2.85 ±0.25 


0.0131 


0.0277 ±0.0031 


1.0 <z< 1.2* 


1.96 ±0.03 


2.55 ±0.20 


0.0162 


0.03 15 ±0.0030 



" Combined COSMOS+COMBO-17 sample 
* COSMOS-only sample 



our power-law fits to small radius would predict that ~ 40% 
of all our galaxies in r < 30 kpc pairs should have r < 15 kpc 
separations; we find that ~ 30% of galaxies in close pairs have 
r < 1 5 kpc separat i ons us ing HST data. Accounting for the fact 
that|Robaina et al. (2009) required additional interaction signa- 
tures such as tidal tails for a pair to be assigned to that cate- 
gorjQ, and that some extremely close pairs will be classified as 
already coalesced merger remnants, we argue that this level of 
agreement is remarkable, and we conclude our extrapolation to 
radii r < 15 kpc is appropriate 

In Table Q] we show the parameters for the 3D correlation 
function that we obtain by deprojecting the 2D 2pcf. The evo- 
lution of the merger fraction is mainly driven by the evolution 
in the number density of galaxies above 5 x 10 10 M Q in physi- 
cal coordinates, together with some evolution in the clustering 
length r . 

In Fig. |2] we show the fraction of massive galaxies found in 
pairs with separations r < 30 kpc as a function of z. To aug- 
ment the data at z ~ 0.1, we calc ulate the SPSS close pair 
fraction using 7 and ro given by iLi et alj d2006l) . We ad- 
just their ro by ~ 5% (an empirical adjustment based on the 
COSMOS+COMBO-17 sample) to account for different lower 
mass limits and binning. We also use the number density of 
galaxies fulfilling our mass criteria, adopting the g-band se- 
lected stellar-mass function in iBellet al.ld2003l) after correcting 
for stellar IMF and H . 

We perform an error-weighted least-squares fit of the form 
F( z ) = /(0) x (1 +z) m to the filled points in Fig. El i.e., the SDSS 
pair fraction, the combined COSMOS+COMBO-17 pair frac- 
tion at 0.2 < z < 0.8 and the COSMOS pair fraction at 0.8 < 
z < 1.2. We find /(0) = 0.0130 ±0.0019 and m= 1.21 ±0.25. 
These associated uncertainties are at the same level as the sys- 
tematic uncertainties in close pair fraction evolution caused by 
uncertainty in the overall stellar M/L scale, and its evolution 
since z = 1 . 

In our correlation function we do not impose a specific mass 
ratio criteria, but given the shape of the mass function above 

2 While most of the projected close pairs at r < 15 kpc are likely to be under- 
going an interaction, the morphological signatures do not need to be necessarily 
evident. Gas fractions and orbital parameters strongly influence the strength of 
the tidal features. Then, we would expect to find a fraction of extremely close 
pairs (respect to the total number of close pairs with r p < 30 kpc) slightly below 
the expected 40%. 

3 Note that we could have seemingly sidestepped this issue by using pairs 
with 15 < r < 30 kpc for the close pair measurement, and adjusting the 
timescale to reflect the amount of time pairs spend separated by 15 < r/kpc < 
30. Nonetheless, in this case, one would still want to confirm that systems that 
reached 15 kpc continue to approach smaller radii, amounting to an exercise 
equivalent to confirming the power law extrapolation. 
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FIG. 1. — Projected 2-point correlation functions for galaxies with masses 
> 5 X 1O 1O M from our combined COSMOS + COMBO-17 sample in dif- 
ferent redshift bins. The best-fit power law is overplotted following the color 
code of the different redshift bins. The power-law parameters are given in 
Table [T] The dot-dashed line represents the b est-fit power law to the SDSS 
sample derived from results by Li et al. ( 2006). The vertical dashed line rep- 
resents the separation limit for close pairs of 30 kpc (physical). The vertical 
dotted line represents the (average of the redshift bins) limit where one can 
no longer reliably resolve individual galaxies in a pair from the ground-based 
photometry. 



M„ > 5 x 10 Mq, we expect most of the mergers to be ma- 
jors; i.e., with mass rations between 1:1 and 1:4. We measure 
the fraction of pairs in projection which fulfill such a mass ra- 
tio criterion, finding a value between 70 and 90 percent in the 
different redshift bins. 

Galaxies from the COSMOS survey represent 70% of our 
sample in the bins ranging from z = 0.2 to z = 0.8, however, the 
addition of COMBO-17 galaxie s helps to decrease the sam- 
ple variance. From the work by Mos ter et alj d2010l) . we esti- 
mate that the sample variance is reduced by a factor of ~ 30% 
by including the three ~ 0.25 sq. deg independent fields from 
COMBO-17. This effect is clearly seen in Fig. [2] Considering 
only galaxies from COSMOS, there is an abrupt transition be- 
tween z = 0.7 and z = 0.9, which is smoothed by the inclusion 
of COMBO-17 galaxies. 

In Fig. [2] we also show the prediction for the close-pair frac- 
tion from the semi-analytic models (SAM) of ISomerville et all 
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FIG. 2. — Fraction of M, > 5 X 10 l0 Mo galaxies in close (3D) pairs 
(r < 30 kpc) as a function of redshift. Diamonds: Pair fraction found using 
only the COSMOS catalog. Black circles: Pair fraction found when adding 
galaxies from the COMBO-17 survey at z < 0.8 (we do not include galaxies 
from the COMBO-17 catalog in the two higher z bins). Star: Pair fraction 
from SDSS. The line shows the best fit to a real space pair fraction evo- 
lution with shape F(z) = /(0) X (1 +z) m , with /(0) = 0.0130 ± 0.0019 and 
m = 1.21 ±0.25 (fit to all black-filled points: star, circles and diamonds). 
Dashed line: Fra ction of galaxies above M * > 5 X 10 10 Mq involved in merg- 
ers, predicted by Somerville et al. 12008) models when assuming a timescale 
of 0.5 Gyrs to transform from merger rate to fraction. 



(2008), updated as described in [Hopkins et ail d2009l) . The 
models s hown here adopt the b est-fit cosmological parameters 
given in Ko matsu et al.l (120091) . We extract the merger rate 
of galaxies where both progenitors have stellar masses above 
5 x 10 I0 M©, normalize by the total number density of galaxies 
above this mass and multiply by two to obtain the fraction of 
galaxies involved in mergers (the quantity that we estimate ob- 
servationally), per Gyr. Then we multiply by the timescale over 
which a merger would be observed as a "close pair", which here 
we assume to be r = .5 Gyr following rPatton & Atfieldl ((2008) 
and iBell et all (2006a) in order to convert this merger rate into 
the equivalent of the observed fraction. With this assumed 
timescale the model predictions are in excellent agreement with 
our observational results. Assuming an uncerta inty of ~ 30% 
in the timescale (matching the numbers from lLotz et al.l 2009 
for close galaxy pairs timescales) would correspondingly shift 
vertically the model line by a similar number. Predictions from 
other SAMs ten d to ag ree reasonably well with those of the 
Som erville et al.l ((2008) models, as shown in Fi g. 1 1 of Jogee 
et al (2009; see also Guo & White 2008 and Hopkins et al.l 
120101) . 

3.1. Comparison with previous works 

Very few studie s of close p air fra ction evolution use mass- 
limited samples. iBelletalJ (l2006al) used COMBO-17 data 
(0.75 sq. deg., the same catalog we use here to complement 
our COSMOS catalog) to find a fraction of galaxies in close 
major pairs (r < 30 kpc) of 2.8% for galaxies more mas- 
sive than 3 x 10 10 M Q at 0.4 < z < 0.8, in excellent agree- 
ment with our result despite the slightly different mass limit. 
They also performed an autocorrelation of all galaxies with 



FIG. 3. — Evolution of the number of galaxies in close pairs. The estimate 
for major pair fraction of galaxies M > 3 X 10 l0 Mp) from Bell et al] <2006al) 
is shown as empty square s. The po int by IXu et'a l. 12004) is shown as the 
grey diamond. Mcintosh et al. (2008) is shown as the grey star. The results of 
Bundv et al. (2009) have been corrected down by 33% in order to match our 
pair definition (see text for details). 



M* > 2.5 x 10 10 M Q finding a pair fraction of ~ 5%. Adopt- 
ing in our analysis instead a mass limit ofM» >2.5x 10 10 M Q 
we recover a pair fraction of ~ 5%. The main driver of this 
higher pair fraction is the fact that ~ 50% of close pairs in a 
sample limited to have M* > 2.5 x 10 10 M Q have mass ratios 
between 1:4 and 1:10; i.e., the major merger fraction is similar, 
but the close pair fraction is boosted by a considerable contri- 
bution from minors. 

Recently, Bun dv et"ai] J2009) performed a similar measure- 
ment by studying a sample of galaxy pairs from the GOODS 
fields (total area ~ 0. 1 sq. deg.). In Fig.[3]we show their results 
for the mass range > 3 x 10 10 M Q after converting their frac- 
tion o f pairs to the fraction of galaxies in pairs. iBundv et al.l 
(2009) study the fraction of major mergers in galaxies with 
M* > 3 x 10 10 M Q , including galaxies 1 .5 mags fainter than that 
mass limit, while by our definition of a close pair, both galax- 
ies need to be above our mass limit. It is natural that they find 
a higher pair fraction, but we quantify this differ ence by per- 
formin g a test in which we try to reproduce the IBundv et al.l 
(2009) methodology. Given the nature of our method it is dif- 
ficult to perform a cross-correlation in order to match their pair 
definition; instead we study the 2D di stribution of galaxies in 
a similar fashion as IBundv et alj d2009l) . That is, we look for 
galaxies with separations r < 30/kpc and matched in redshift 
space and subtract random counts. We perform such a test in 
the redshift range 0.2 < z < 0.6, where we are complete for 
stellar masses of ~ 1 - 2 x 10 10 M Q , with both our pair defini- 
tion (both galaxies above 5 x 10 10 M Q ) and also keeping such 
a mass limit for the primary but allowing the secondary to be 
a factor 4 less massive. In the first case we find, as expected, 
a result compatible with the one found with the 2pcf, and in 
the second we find that the pair fraction is larger by a factor of 
~ 1.5. Extrapolating th is ratio to the higher redshift bins, we 
accordingly re-scale the IBundv et alj (120091) numbers by mul- 
tiplying their results by 0.67 in order to match our pair defini- 
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tion. The difference is reduced to a factor of less than 2 in the 
lower redshift bins, while our re sults basi c ally a gree at z ~ 1 . 
We also show the estimate from IXu et al.l d2004l) . who used a 
combined sample from the 2MASS and 2dFGRS surveys, af- 
ter c onverting their r e sults t o our IMF and Hq, and the result 
from McIntoshetaD d2008), who measured the merger frac- 
tion in galaxy groups based in projected pairs with additional 
signs of gravitational interaction. Before performing the com- 
pariso n, we further corre c t the p air fractions found bv lXu et ail 
(120041) and iBundv et all d2009h down by 30% to account for 
pairs in their analyses that are genuinely associated with each 
other (so are not random projections), have projected separa- 
tions of < 30kpc but are separated by more than 30kpc in real 
space (i. e., galaxies in g roups that are projected along the line 

of sight; EnetaDllOOMQ- 

We can not compare our measurements with morphological 
studies of merger fractions because of uncertainties with the 
nature of the progenitors ( differentiation between major and 
minor mergers can be extremely challenging, especially in the 
case in which one of the progenitors is very gas-rich), and selec- 
tion effects related to orbital parameters, galaxy structure and 
gas fractions, with the latter introd ucing a bias in th e recovery 
of interacting systems at higher z (lLotz et al.ll2009b . A more 
fundamental difference with respect to morphologic al studies, 
both those performed with visual clas sifications (jJogee et al.l 
20091) and with automated classifiers ( Conselice et alj 120091: 
Lotz et al. I2008I) . is the fact that both merger remnants and 
closely interacting pair s can show si gnatures of gravitational 
interactions (see [Robaina et al. 2009, for a discussion on the 
relative numbers of merger remnants and close pair interactions 
with morphological disturbances). These two groups are related 
to different time scales, and the strength and abundance of tidal 
features is also different (often degenerate with other parame- 
ters such as gas fraction). Studies on close pair statistics aim to 
detect the galaxy interaction in an early stage, where both pro- 
genitors can be studied separately, while morphological studies 
include systems both before and after the coalescence. For this 
reason, a direct, quantitative comparison is impossible. It is 
also hazardous to compare with cl ose pair measu rements based 
on luminosity-selected samples. iLin et alj d2008l) used a sam- 
ple of galaxies with B-band magnitudes (corrected for passive 
luminosity evolution) -21 < Mb < -19. For red galaxies, this 
is roughly compatible with our mass-selected sample, but it 
includes many low-mass blue galaxies, which makes a com- 
parison impossible because their clustering properties are very 
different. As an extreme example, in which one close pair is 
formed by a very red Mg = -21 galaxy and a very blue one with 
Mb = -19, the stellar-mass ratio in this hypothetical interaction 
could be as high as 40. This pair definition is clearly not com- 
patible with one in which both galaxies are selected by stellar 
mass. Furth ermore, as merging can enhance the s tar forma- 
tion activity dBarton et al.ll2000URobaina et al.ll2009h . selecting 
galaxies in rest-frame blue bands is biased in favor of mergers, 
as such a selection would recover merging systems with lower 
mass given the decreased M/L ratio induced by the interaction. 
All of these effects will be more pronounced in gas-rich galax- 
ies, and gas fractions were likely higher in the past, leading to 
a redshift dependent bias . 

iKartaltepe et al.l d2007l) derived the evolution of the pair frac- 
tion from a luminosity-selected sample from a similar dataset 

4 This effect is present also in pair fractions determined from spectroscopic 
redshifts as long as a deprojection to the 3D space is not performed. 



as analyzed in this paper. Their very different result (they find 
m = 3.1 ± 0.1) is caused at least in part by the difference be- 
tween mass- and luminosity-selected samples, as described 
above, as well as by the fact that they do not correct for pas- 
sive luminosity evolution. In addition, they identify pairs in 
both the ground-based and HST/ACS-based catalogues. Be- 
cause their sample is selected by ground-based luminosity, very 
close pairs that are only resolved by ACS can be as bright as a 
single galaxy in a more widely separated pair that is resolved 
in the ground-based imaging. This artificially raises that close 
pair fraction, especially at high redshift. 

Summarizing, we choose not to compare with merger frac- 
tions obtained by morphological indicators, nor with those ob- 
tained by close pair counts in luminosity-selected samples; 
these provide measurements of different merger-related quan- 
tities and can not, in general, be compared to our measurement. 

3.2. The impact of galaxy merging on the creation of 

M* > 10 n M Q galaxies. 

Our measured fraction of > 5 x 1O 1O M galaxies in 
close pairs as a function of redshift constrains the impact that 
merging-induced mass assembly has on the creation rate of 
> lO n M systems. 

As only ~ 6% of our galaxies with a projected companion 
at r < 30 kpc has a second companion at such separation, we 
assume that the number of close pairs is simply one half of the 
number of galaxies in close pairs, neglecting the very minor 
influence of such a small number of galaxies in triplets. The 
fraction of newly created M* > 10 u Mq galaxies due to merg- 
ing is f rem = N cp /N\i, where N cp is the number of close pairs of 
galaxies above 5 x 1O 1O M each (regardless of their optical col- 
ors) and N\ \ is the total number of galaxies with stellar masses 
in excess of lQ ll M ( r l . 

Following lPatton & Atfieldl d2008l) or lBell et all d2006al) . the 
merger timescale for galaxy pairs at this separation is approxi- 
mately t = 0.5 Gyrs, so we compute the creation rate of newly 
assembled galaxies, R rem , as R rem = f rem /T. We integrate the 
merger rate over cosmic time finding that, on average, present 
day galaxies with stellar masses larger than 10 n M Q have un- 
dergone 0.5 (0.7) mergers since z = 0.6 (1.2) from interactions 
between galaxies more massive than 5 x 10 10 M Q . 

3.3. The merger-driven evolution of the red sequence. 

In the previous sections we have studied the impact of merg- 
ing of the evolution of the mass function of galaxies, without 
restriction to the optical colors of the progenitors or descen- 
dants. Now we want to address the question whether our ob- 
served merger rate evolution can explain the observed number 
density evolution in the massive end (M„ > 10 u Mq) of the red 
sequence. 

First, we convert the evolution of the B-band luminosity 
function (LF) of red galaxies as measured by Brown et al. 
(2007) to the evolution in the number density of galaxies more 
massive than 10 u Mq. For this conversion we use the typi- 
cal stellar mas s-to-light ratio of nearby red galaxies ( M/Lb = 
3.4 ±0.6, e.g jBell et all2003llKauffmann et alj2003l) and ap- 
ply a correction to account for its evolution with redshift, de- 
rived from the evolution of the fundamental plane (FP) zero 
point (Alog(Af /Lb) = 0.555 ± 042 between z = 1 and z = 0; 
Ivan D okkum & van der Marell 120071) . We repeat the process 
10000 times allowing the luminosity function parameters, as 
well as the M /Lb constraints, to vary randomly within their un- 
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certainties. We estimate the final uncertainty by estimating the 
typical dispersion of those 10000 realizations. The results from 
this exercise are shown in Fig. [4] We find an evolution in the 
number density of the massive end of the red sequence of a fac- 
tor ~ 4 between z = 1 and z = 0. 

Our inference of a factor ~ 4 evolution in the number density 
of red sequence galaxies with M* > 1 1 1 Mq is somewhat larger 
than a number of recent estimates. The main reason for this is 
our ad option of the recent results by van Dok kum & van der Marell 
d2007l) for the evolution of M/L determined through FP evo- 
lution, who find relatively rapid evolution of M/L leading to 
more rapid inferred fading of the stellar popul ation than as- 
sumed by many previous works. For example, ICimatti et alj 
(2006) found a result compatible with a constant number den- 
sity of M* > 10 n M Q galaxies in the interval < z < 0.8. In 
order to obtain that result they applied a passive luminosity 
correction of Alog(M/Ls) = 0. 46 ± 0.04 to a sample drawn 
from COMBO-17 and DEEP2 (iDavis et al J 120031: iFaber etafl 
2007 ;). This luminosity correction was obtained from the works 
bv lvan Dokkum & Stanford! d2003l) and lTreu et ai]d2005l) and is 
lower than the one found a few years later bv lvan Dokkum & van der 
d2007h with a larger sample. This apparently small differ- 
ence in Alog(M/Lfi) has strong implications when studying the 
number density evolution of massive red galaxies at the level 
of ~ 1.5-2 between z = 1 and z = 0. Also, the number of 
galaxies in their COMBO-17 or DEEP2 samples is relatively 
small when compared to the ~ 7 sq. deg. of the NOA O Deep 
Wide-Field Survey (NDWFS, llannuzi & Devi 1 19991) used by 
iBrown et al! (120071) . ICool et all d2008l) used asample of lumi- 
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FIG. 4. — Number density evolution of red galaxies with M„ > \0 11 Mq. 
Filled point s with error bars represent the number density observed by 
when translated to stellar mass (see text for details). The 



nous red galaxies (LRG) from SDSS and spectroscopic follow- 
up from MMT to address a similar question. They found no 
evolution in the number density of LRGs since z = 0.9. We be- 
lieve that the main driver of the difference between Cool et al.'s 
result and ours is, again, the correction for M/L ratio evolution. 
They used stellar populatio n models from iBruzual & Chariot! 
(2003) and Mara stonl d2005l) . which gives slower evolution than 
FP zeropoint evolutio n constraints. 

In a recent work, [Matsuoka & Kawaral d2010h studied the 
number density evolution of massive galaxies since z ~ 1 from 
a sample drawn from UKIRT Infrared Digital Sky Survey and 
SDSS. Instead of applying an average M/L correction, they cal- 
culated stellar masses for individual galaxies at all redshifts, 
finding a factor ~ 4 evolution at M* > 1O U M0 (when adding - 
up their two mass bins). Unlike our present work, they did 
not differenciate between red and blue galaxies, but given the 
predominance of red galaxies at such high stellar masses (they 
report a blue fraction < 30% at all redshifts in the same mass 
range) we believe that our factor ~ 4 evolution agrees remark- 
ably well with their result. 

We now address the question of whether galaxy mergers be- 
tween massive galaxies can drive the observed number density 
evolution of massive red galaxies. For that purpose, we use the 
result for the merger rate found in Section 13.21 and work un- 
der the assumption that mergers between massive galaxies will 
produce remnants with red optical colors. Dry (gas-free) merg- 
ers are likely to play an importa nt role in the build-up of the 
massiv e end of the red sequence (Ivan Dokkum 2005; Bell et al. 
2006b), and the descendant systems of such interactions will 
clearly be red systems as well. Moreover, mergers involving 
massive blue, disk-like galaxies are also expected to produce 
red ellipticals as remnants. Major interactions, independently 
of the color of the progenitors, are expected to perturb the or- 



et all 120071) 

solid Tine shows the expected growth implied by our measurement of the 
merger rate. 



bits of the stars in the galaxies and form elliptical systems, 
but the observed fraction of spheroidal galaxies with blue opti- 
cal colors is a negligible fraction of the total at masses above 
10"M^ (see [Schawins kilitail 120091: Ikannappan et all 12009: 



iHuertas-Companv et al.ll20Tol) . It is likely that even those few 
systems will passively evolve to the red sequence in < 1 Gyr, 
so, our assumption that merger remnants with masses above 
such stellar mass will become red sequence galaxies is justi- 
fied. 

We use the number density of massive systems at z ~ 0.9 
that we have obtained as the starting point and calculate the 
growth in the number of massive red galaxies implied by our 
measurement of the merger rate. We show the result of this ex- 
ercise as the solid line in Fig. [4] We stress that the observed 
evolution (filled circles) and the one implied by our measure- 
ment (solid line) are completely independent except for the fact 
that we use the observed density at z ~ 0.9 to anchor the evolu- 
tion predicted by our close pair fractions. We find that mergers 
of massive galaxies can explain the evolution in the observed 
number density of massive red galaxies since z = 1 . In different 
words, the majority of galaxies observed today at the massive 
end of the red sequence have been assembled in the last 8 Gyrs 
through mergers of massive galaxies. 

W e have used r = 0.5 Gy rs, but using the t ~ 1 Gyr timescale 
from Kitzbich ler & White! d2007l) produces a somewhat slower 
evolution that is still compatible within the error bars. 

There are two caveats we would like to mention. Firstly, 
given the nature of our method, some of the progenitor galaxies 
have masses above 1O H M0, so strictly speaking they are not 
newly formed massive galaxies. Second, because we adopt a 
lower mass limit of 5 x 10 10 M Q , we underestimate the num- 
ber of major mergers that could lead to the formation of a 
massive galaxy. For example, a pair with individual masses 
M* = 6 x 1O IO M and M„=4x 10 10 M Q would not make it into 
our pair sample but would produce a major merger remnant of 
1O U M0. Assuming that the merging population has a compo- 
sition similar to our overall sample, we estimate that this latter 
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lower mass limit issue has an impact a factor ~ 2 larger on the 
creation of > 10 n M Q galaxies than the overestimate caused by 
double counting already massive red galaxies (i.e., for every 10 
massive galaxies we are counting as produced by mergers, there 
would be 1-2 which were already above 1O U M and ~ 3 that 
we are not accounting for because one of the progenitors has a 
mass below our mass limit of 5 x 10 10 M Q ). This means that our 
results are uncertain at the level of < 20% when counting newly 
formed massive red galaxies. We do not apply this correction to 
our results but would like to stress that, if anything, we slightly 
underestimate the real creation rate of massive galaxies induced 
b y merging. 

iHopkins et all d2008l) found by using HOD models a merger- 
driven evolution of a factor 3-4 in the integrated stellar mass 
density in red elliptical galaxies between z = 1 and z = 0. Given 
that we focus on the massive end of the distribution, a quanti- 
tative comparison is impossible, but we argue that our results 
are in qualitat ive good ag r eemen t with the trend seen in that 
work and in Hopki ns et alj (|2006), where they found some evo- 
lution in M* > 10 u Mq red galaxies in the same redshift range 
by evolving the mass function of red galaxies in a manner com- 
p atible with th e predi cted merger rates. 

IWhite et alJ (2007) studied the clustering of luminous red 
galaxies in the NDWFS in conjunction with HOD models, find- 
ing that ~ 1 /3 of red luminous satellites present at z = 0.9 have 
dissapeared by z = 0.5. The most likely explanation to this phe- 
nomenon is that merging has taken place, decreasing the num- 
ber of such galaxies, in good qualitative agreement with our 
result. 

Also, in a recent work, Ivan der Wei et al. I (120091) studied the 
shape (axis ratio) distribution of a large sample of quiescent 
galaxies drawn from the SDSS. Quiescent galaxies with promi- 
nent disks are exceedingly rare at stellar masses above lO n M , 
providing a strong indication that most, if not all, massive qui- 
escent galaxies form through major mergers. Our results are 
consistent with this entirely indepent result, and, moreover, sug- 
gest that most of this merger activity occurred over the past 
~ 8 Gyrs, since z = 1 . 

4. CONCLUSIONS 

We have studied the impact of galaxy mergers on the evolu- 
tion of massive galaxies by using 2-point correlation functions 
to measure the fraction of galaxies in close pairs from a sam- 
ple of ~ 18000 galaxies more massive than 5 x 10 10 M Q drawn 
from the COSMOS and COMBO-17 surveys and a pair frac- 
tion estimate from SDSS. We have also used constraints from 
the observed evolution of the fundamental plane zeropoint to 
calculate the number density evo lution of M* > 10 n M Q red 
galaxies from lBrown et al. (2007) LFs. Our main findings are: 

1- . The fraction of galaxies more massive than 5 x 10 10 M o in 
close pairs evolves as F(z) = (0.0130±0.0019) x (l+z) 121±a25 
in the redshift range < z < 1.2. Assuming a merging timescale 
of t = 0.5 Gyr this implies that galaxies more massive than 
10 n M Q have undergone, on average, 0.5 (0.7) major merg- 
ers since z = 0.6 (1.2) involving progenitors more massive than 
5 x 1O 1O M . 

2- . The number density of M* > 10 h Mq galaxies on the 
red sequence increases by a factor of ~ 4 between z = 1 and 
z = 0. This result depends strongly on the assumed evolu- 
tion of the M/L ratio, which we have assumed here to be 
A log(M/L B ) = 0.555 ±0.042 between z= 1 and z = following 



recent estimates by van Do kkum & van der Marell d2007l) . 



3-. The evolution implied by our measured merger rate is 
sufficient to explain this observed number density evolution of 
massive red galaxies since z = 1; major mergers between mas- 
sive (M* > 5 x 10 10 M Q ) galaxies are the main driver of the red 
sequence growth at its massive end. 
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